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Abstract. We present first elements of kinetic theory appropriate to the inhomogeneous phase of the 
Hamiltonian Mean Field (HMF) model. In particular, we investigate the case of strongly inhomogeneous 
distributions for T — > and exhibit curious behaviour of the force auto-correlation function and friction 
coefficient. The temporal correlation function of the force has an oscillatory behaviour which averages to 
zero over a period. By contrast, the effects of friction accumulate with time and the friction coefficient 
does not satisfy the Einstein relation. On the contrary, it presents the peculiarity to increase linearly with 
time. Motivated by this result, we provide analytical solutions of a simplified kinetic equation with a time 
dependent friction. Analogies with self-gravitating systems and other systems with long-range interactions 
are also mentioned. 

PACS. 05.20.-y Classical statistical mechanics - 05.45.-a Nonlinear dynamics 



! 1 Introduction 

■ Recently, there was a renewed interest for the statisti- 
' cal mechanics of systems with long-range interactions 0. 
, This concerns self-gravitating systems (galaxies) in as- 

■ trophysics |2I3| . large-scale coherent structures (jets and 
' vortices) in geophysical flows bacterial populations 
. (chemotaxis) in biology |4I5| . clusters in the Hamiltonian 
' Mean Field (HMF) and Brownian Mean Field (BMF) mod- 
, els |til7| . galactic bars !Sj, neutral and non- neutral plasmas, 

' dislocation dynamics, planetary formation, cosmology etc. 
' The dynamical evolution of such systems presents a lot of 
I peculiarities (9j. For Hamiltonian systems with long-range 
' interactions, the collisional relaxation time diverges with 
, the number N of particles so that the system experiences 
two successive types of relaxation: a collisionless relax- 
ation on a short timescale of the order of a few dynamical 
. times to (called violent relaxation in astrophysics) and 
' a collisional relaxation on a long time scale of the order 
NHd with S > 1. The first regime leads to the forma- 
tion of a metaequilibrium state, or quasi-stationary state 
(QSS), which is a stable stationary solution of the Vlasov 
equation that is not necessarily of the Boltzmann form 
jlO 11 12 13 14 15 16 . The second regime leads in general 
to the ordinary statistical equilibrium state described by 
the Boltzmann distribution. In the case of self-gravitating 
systems, there may not exist statistical equilibrium and 
the system can evaporate or collapse (gravothermal catas- 
trophe) [2]. Between the phase of violent relaxation and 
the late collisional evolution (equilibrium or collapse) the 
system passes by a succession of quasi-stationary states 



which are quasi-stationary solutions of the Vlasov equa- 
tion slowly evolving with time under the effect of en- 
counters (flnite N effects). In astrophysics, this phase is 
described by the orbit-averaged-Fokker-Planck equation 

ri8-i4 - 

The developement of a kinetic theory for the collisional 
evolution of systems with long-range interactions is com- 
plicated for different reasons. First of all, the temporal 
correlation function of the force may not decay sufffciently 
rapidly to vindicate the Markovian approximation that is 
used in many kinetic theories. For example, in stellar dy- 
namics, the temporal correlation function decreases like 
t~^ leading to a logarithmic divergence of the diffusion 
coefficient ^Hl- In the case of the HMF model, the cor- 
relation function decreases exponentially rapidly but the 
correlation time diverges close to the critical point T Tc 
|20l7j . This is a general feature of long-range attractive 
potentials of interaction On the other hand, systems 
with long-range interactions are usually spatially inhomo- 
geneous and non-local effects strongly complicate the ki- 
netic theory. In the case of self-gravitating systems, one 
usually avoids the problem by making a local approxima- 
tion and developing the kinetic theory as if the system 
were homogeneous ^ISj. This is partly justified by the fact 
that the fiuctuations of the gravitational force are dom- 
inated by the contribution of the nearest neighbor (the 
distribution of the force is a particular Levy law called 
the Holtzmark distribution) On the other hand, for 
the HMF model, the kinetic theory has been developed 
only in the case T > Tc where the system is in a sta- 
ble homogeneous phase |22l7l23j . In these situations the 
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meanfield force vanishes and, to a first approximation, the 
particles follow linear trajectories with constant velocity. 

One goal of this paper is to present elements of kinetic 
theories valid for the inhomogeneous phase of the HMF 
model (T < Tc) and show that the situation becomes sen- 
sibly different from what we are used to in the case of 
homogeneous systems. In particular, we shall investigate 
the case of strongly inhomogeneous systems for T ^ 0, 
where the particles cluster around 9 = 0. In that case, 
their mean motion is that of a harmonic oscillator and 
it is possible to calculate analytically the auto-correlation 
function of the force and the friction. We find that these 
quantities present a curious behaviour. The temporal cor- 
relation function of the force has an oscillatory evolution 
which averages to zero over a period. By contrast, the 
effects of friction accumulate with time and the friction 
coefficient does not satisfy the Einstein relation. On the 
contrary, it presents the peculiarity to increase linearly 
with time. These curious behaviours were previously noted 
by Kandrup '1 1 in the case of self-gravitating systems but 
the consideration of the HMF model allows to obtain more 
explicit results (devoid of any gravitational divergences) 
and provides a simple framework where these effects can 
be studied in detail. 

The paper is organized as follows. In Sec. |21 we re- 
call basic results concerning the structure of the statistical 
equilibrium states of the HMF model. In Sec. 01 we develop 
a kinetic theory of the HMF model valid for both the ho- 
mogeneous and the inhomogeneous phase. We present a 
generalized Landau equation describing the evolution of 
the distribution function of the system as a whole and a 
generalized Fokker-Planck equation describing the evolu- 
tion of the distribution function of a test particle (or an 
ensemble of non-interacting test particles) in a thermal 
bath of field particles at statistical equilibrium. In Sec. 01 
we calculate the temporal auto-correlation function of the 
force acting on the test particle. We show that it presents 
an oscillatory behavior which averages to zero over a pe- 
riod. We also discuss the expression of the diffusion coef- 
ficient and its relation to the Kubo formula. In Sec. El we 
calculate the friction force acting on the test particle. We 
show that the frictional effects are cumulative and lead to 
a linear divergence of the friction coefficient for t +oo. 
We also discuss the break-up of the Einstein relation for a 
strongly inhomogeneous system. In Sec. El we provide the 
analytical solution of a simplified kinetic equation with a 
time-dependent friction related to our study. 



2 Statistical equilibrium states of the HMF 
model 

The HMF model consists of N particles (of unit mass) 
moving on a ring and interacting via a cosinusoidal poten- 
tial. The phase space coordinates {Oi,Vi) of the particles 
satisfy the Hamiltonian equations of motion 



dt 



dH 

dvi ' 



'dt 



dH 
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^ 1 



1=1 



k 

An 



^cos( 



(1) 



The HMF model was introduced by several groups (see a 
short historic in 7 ) either as a simple model with long- 
range interactions mimicking gravitational dynamics 25 261^ 
or as a simplified model for the formation of bars in disk- 
shape galaxies jSl. The thermodynamic limit corresponds 
to iV ^ -\-oo in such a way that the rescaled temperature 
f] = kM/inT and rescaled energy e = 87r£'/fcM^ remain 
of order unity (this can be conveniently accomplished by 
taking the coupling constant k ^ 1/N right from the be- 
gining, in which case T ~ 1 and E ^ N). In this proper 
thermodynamic limit, the mean- field approximation be- 
comes exact for N —^ +oo, except near the critical point 
At statistical equilibrium (see, e.g., |7]), the distribu- 
tion function can be written f{9, v) = p(9)f{v) with 



1/2 



p{0) 



M 



2nIo{m 



-0B cose 



where we have adopted the normalization J^^ p{9)d9 
and f{v)dv = 1 (note that here M = Nm — N). 
Equation ^ is the Boltzmann distribution in a potential 
<P = B cos 9. The meanfield force experienced by a particle 
is 



(2) 



(3) 



M 



where 



{F){9) = -'P\9) = S sin 61, 



B = -— / p{9' ) cos 9' d9'. 
27r Jo 



(4) 



(5) 



To obtain Eqs. Q-©, we have assumed (without loss of 
generality) that the equilibrium distribution of the system 
is symmetric with respect to the a;-axis. The quantity B 
is similar to the magnetization in spin systems. It is de- 
termined as a function of the temperature (see e.g. |7]) by 
the implicit equation 



B 



kMh{j3B) 
27r loim' 



(6) 



The energy is given by = NT 12 - nB'^/k. For T > 
Tc = kM/An 01 E > fcM^/STr, the only solution to 

the above equation is i? = leading to an homogeneous 
distribution of particles. For T < Tc or E < Ec, the homo- 
geneous phase becomes unstable and a (stable) clustered 
phase appears with B 0. This corresponds to a second 
order phase transition (see e.g. \7\). 



3 The inhomogeneous kinetic equation 

We shall discuss here the kinetic theory of the HMF model 
by using general results coming from the projection oper- 
ator formalism. This formalism starts from the Liouville 
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equation for the A^-body distribution function P^dOi, Vi}, t) 
and derives an exact kinetic equation for the one-body dis- 
tribution function f{0,v,t) — NPi {6, v,t) by using projec- 
tion technics. This equation is then simpUfied by making 
some approximations on the correlation function of the 
field particles. This formalism introduced by Willis & Pi- 
card 28 is quite general and leads to a form of general- 
ized Landau equation 9, . It was applied by Kandrup (29j 
in the case of stellar systems, by Chavanis |30l3j for two- 
dimensional point vortices and in 7 for the HMF model. 
This formalism is also very close to the linear response 
theory developed in |24l31j where the friction term (or 
the drift term in the case of point vortices) is calculated 
directly from the perturbation of the A^-body distribution 
function of the bath caused by the interaction with the 
test particle. One interest of this formalism is that it re- 
mains valid in the case of non-Markovian and spatially 
inhomogeneous systems while other formalisms developed 
in connection with plasma physics (binary collision mod- 
els, BBGKY hierarchy, diagrammatic expansions, quasi- 
linear theory,...) usually consider Markovian and homoge- 
neous systems and work in Fourier space. By contrast, the 
projection operator formalism remains in physical space 
which enlightens the basic physics. We shall not repeat the 
intermediate steps of the formalism which can be found in 
128129130] . This formalism leads to a general kinetic equa- 
tion of the form 



dt d9 ^ ' dv dv 



J dr J d0idviT{l 0,t) 



xG{t,t-T) 



HI 



0,t-T)—+T{Q 
ov 



xf{e,v,t^T)f{0,,v,,t-T), (7) 



where G(i, t — t) is the Green function associated to the 
averaged Liouville operator constructed with the mean- 
field force {F){d, t) and JF(1 — > 0, t) is the fluctuating force 
created by particle 1 on particle (see, e.g., j22. for more 
details). Under this form, we clearly see the terms of diffu- 
sion and friction (first and second terms in the r.h.s. of Eq. 
JTJ) and their connection to a generalized form of Kubo 
formula (the time integral of the force auto-correlation 
function). These points will be developed in the sequel. 
The ratio of the right hand side (collision term) on the 
left hand side (meanfield advective term) is of order 1 /N 
in the proper thermodynamic limit |7I9| recalled in Sec. 
12 Therefore, for N — > +oo, Eq. {T)) reduces to the Vlasov 
equation. The collision term takes into account finite N 
effects and can be viewed as the first order correction to 
the Vlasov limit in an expansion in A^"-'^ (see, e.g., |2| p. 
260). Equation ((JJ) can be viewed as a generalization of the 
Landau equation (initially introduced in plasma physics) 
to which it reduces ^ if the system is homogeneous and 
Markovian One drawback, however, of the projection 
operator formalism (or more precisely of the approxima- 
tions leading to Eq. (01) is that it ignores collective ef- 

^ Note, however, that for one-dimensional systems such as 
the HMF model, the Landau collision term cancels out |22I7I9| : 
see also last paragraph of |32| . 



fects which are important especially close to the critical 
point. Such collective effects can be taken into account by 
using the Lenard-Balescu equation in the case of homo- 
geneous systems |22l7l9j . We shall not discuss these col- 
lective effects here and shall remain close to the situation 
considered by Kandrup in the astrophysical setting 
by adapting and expliciting the calculations in the case of 
the HMF model. 

Equation ((TJ is an integrodifferential equation (with 
respect to the variables 9i, vi) describing the evolution 
of the system as a whole. We shall consider here a sim- 
pler problem, namely the evolution of a test particle in 
a bath of field particles with prescribed static distribution 
f{Oi, vi) which is a stable stationary solution of the Vlasov 
equation. By adapting the projection operator formalism 
to this situation where f{9i, vi) is fixed, we find that the 
time evolution of the density probability P{9, v, t) of find- 
ing the test particle in [9, v) at time t is governed by the 
equation 



dP dP dP 9 /■* f 



t) 



HI 



o,t-T)— + Ho 

ov 



xP{9{t -T),v{t~T),t~ T)f{9l {t-T),V,{t~T)). (8) 

In this paper, we shall consider the evolution of a test 
particle in a thermal bath of field particles at statistical 
equilibrium described by the distribution ©-JSl- that 
case, we obtain a kinetic equation of the form 



dP dP 



dt 



89 



{F) 



dP _ d 
dv dv 
d 



dr J dOidviTil ~* 0, i) 



x<^ J^(l ^ 0,f-r)— l,t-T)/3vi(t-T) 
I ov 

xP{9{t ~ T),v{t ~T),t~ T)p{0,)f{v^). (9) 

The fiuctuating force can be written T{1 ^ 0,t) = F{1 — > 
0, t) - (F(l 0, t)) where 



F(1^0,i) 



2tt 



smi9{t)~9i{t)), 



(10) 



is the exact value of the force created by particle 1 on 
particle and 



(F(1^0,i)) 



1 

N 



F{1^0,t)p{9i)f{vi)d9,dvi, 



(11) 

is its mean-field value. Equation can be seen as a sort 
of generalized Fokker-Planck equation. However, the dy- 
namics is generally non-Markovian (see below) so that Eq. 
is not a Fokker-Planck equation. 

The time integral in Eq. ij^ must be performed by 
moving the particles with the meanfield force between 
t and t — T (see, e.g., jSH])- Accordingly, the quantities 
9{t~T) and v{t—T) are solutions of the equation of motion 



f9 
dt^ 



B sin 9. 



(12) 
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This is the equation of a nonlinear oscillator. The general 
solution is given by 



e{t) 



y/vQ + 2i?(cos 9o ~ cos ( 



±t, 



vit) 



± 



2S(cos6'o - cos9{t)). 



(13) 



(14) 



In our previous study we have considered the case 
T > Tc where the distribution of the bath is homoge- 
neous and the particles follow linear trajectory with con- 
stant velocity in a first approximation. In that case, Eq. 
© reduces to a Fokker-Planck equation of the Kramers 
form with a friction coefficient given by the Einstein rela- 
tion. This kinetic equation converges for t —^ +oo to the 
Maxwellian distribution of the bath |7]. However, the dif- 
fusion coefficient decreases rapidly with the velocity (like 
the Gaussian distribution of the bath) and this leads to 
anomalous diffusion ^22, and to a slow progression of the 
front in the high velocity tail [221 ■ Here, we shall consider 
the case T ^ where the bath distribution is strongly 
inhomogeneous. Assuming B — — cj^ < (the case B > 
is symmetric), the particles cluster around 9 = and form 
a Dirac peak for T = 0. For T ^ 0, they will remain lo- 
calized around 9 = 0. Thus, we can expand Eq. H12(l to 
first order in 9. It becomes the equation of a harmonic 
oscillator 



and the equations of motion are explicitly given by 
9{t) = — sm{ujt) + 6*0 cos(wt), 

UJ 



v{t) = vq cos{ujt) — 9quj sin(wi). 
For T ^ Tc, the parameter B is given by 



B 



Bn 



= ± 1- 



T 



(15) 



(16) 



(17) 



(18) 



where B^ax = ^Tc = kM/2'K |7). Therefore, to leading 
order, we find that the pulsation of the particle trajectory 



IS 



kM 
"2^' 



(19) 



We shall not attempt in this paper to study the full 
kinetic equation l^. We shall limit ourselves to analyti- 
cally study the behaviour of certain quantities that enter 
in this equation. To emphasize the physical structure of 
the kinetic equation we introduce two functions 

C(t, t)= j d9^dviT{l ^ 0, t)T{\ ^ 0, t - T)p{9i)f{vi), 

(20) 



^{t, ^) = /3 y d9idviT{l 0, t) 

^ 0, i - T)vi{t - T)p{9i)f{vi), (21) 
in terms of which Eq. lO can be rewritten 

-<F{t, T)}p{9{t -T),v{t~T),t-T). (22) 



The function C(t,T) which appears in the diffusion term 
represents the force auto-correlation function. It will be 
specifically studied in Sec.0] On the other hand, the func- 
tion if'(t, r) which appears in the drift term is connected 
to the friction force which will be specifically studied in 
Sec. |5l (see in particular the quantity /). Then, the diffi- 
culty with Eq. (|22|l lies in that fact that it is non Marko- 
vian: we need to know the value of P(9,v,t) at previous 
times to pass from t to t -I- dt. In Sec. EJ in order to get 
some insight into the problem, we shall consider a simpler 
kinetic equation where we neglect non-Markovian terms 
altogether but keep the time dependence of the diffusion 
and friction coefficients. 



4 The force auto-correlation function 

4.1 The temporal correlation function 

One quantity of great interest in the kinetic theory is the 
force auto-correlation function. Indeed, in ordinary cir- 
cumstances, the diffusion coefficient in the Fokker-Planck 
equation is expressed as the time integral of the auto- 
correlation function through the Kubo formula. The tem- 
poral auto-correlation of the fluctuating force can be de- 
composed into 

(mm) = (mm) - {m){m), m 

where J^(i) = F{t) ~ {F{t)) is the total fluctuating force 
acting on the test particle at time t. Using F{t) = P{i ~^ 
0, t), we get 

(F(O)F(t)) - ^(F(z ^ 0, 0)F{j ^ 0, t)) 

id 

= ^(F(i->0,0)F(i->0,i)) 

(*^0,0)F(j->0,t)). (24) 



Since the TV-body distribution of the bath is a product 
of N one-body distributions (see, e.g., [7]) and since the 
particles are identical, we obtain 

(P(O)F(t)) = N{F{1 -> 0,0)F(1 ^ OA)) 
+N{N - 1)(F(1 ^ 0,0))(F(1 ^ 0,t)). (25) 
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Accounting that {F{t)) = N{F{1 0,t)), we get 

(HOmt)) = N{F{1 ^ 0, 0)^^(1 ^ 0, t)) 
-N{F{l->0,0)){F{1^0,t)), 

which can be written 

(^(O)^(t)) = iV(^(l ^ 0, 0)^(1 ^ 0, i)), 



In the preceding expansions, we have implicitly assumed 
that the coordinates of the field particles and of the test 
particle scale as w ~ •\/2//3 and 6 ~ (l/w)-\/2//3. There- 
(26) fore, our asymptotic expansion is valid to order T for 

Now, the correlation function of the fluctuating force 



(27) 



is given by 



where ^(1 - 0, i) = ^^(1 -> 0,t) - (^^(1 ^ 0,t)) is the (mm) = {^(O)F(i)} - iV(F(l ^ 0,0)) (^^(1 ^ 0,^)), 
fluctuating force produced by particle 1 on the test parti- (36) 
cle. Explicitly, 



where 



(28) 



27riV 



sm{eit) ~ di{t))p{9i)fivi)d0idvi. 

(37) 



We note that this combination of terms enters in the dif- 
fusion term in Eq. JH)!; this is precisely the function C{t, r) 
defined in Eq. (|20|l where we have taken the origin of times 
at t = 0. Let us first compute the quantity 



Using the same approximations as before, we obtain 



{F{1 0,t)) = — — [-sin(wt) +ecoa{Ljt) 



{F{0)F{t)} = N(F{1 ^ 0,0)F(1 ^ 0,t)) 
Explicitly, we have 

{F{0)Fit)} ^ j F{1^ 0,0)F(1 ^ 0,i) 

xp{e^)f{vi)de^dvi. 

Using Eq. lfTU|l . we get 



(29) Combining the previous results, we get 

fc2 M 



(mm) 



47r2 /3w2 

Using Eq. (|19|l . we finally obtain 



cos(ti;t). 



(30) 



(•^(0)m = ;^cos(c^0 + O(T2)- 
Z7rp 



(38) 



(39) 



(40) 



{mm} 



47r2 



0i)sin(0(t)-0i(t)) 
xp{ei)f{vi)deidvi, (31) 



We note that, to order T, the correlation function of the 
fiuctuating force depends only on the ellapsed time t and 
not on the initial instant t = 0. We also note that the 
correlation function is periodic with the same pulsation co 



where 9,(1) denotes the position at time t of the i-th par- the particle trajectory and that it averages to zero over 
tide located at 9i at t = 0. Now, using the equation of ^ period, 
motion ((TB|l . we obtain 



9{t) - 9i{t) = - sui{ujt) + (j>cos{ujt), 



(32) 



where (t> = 9-9i&ndu = v-vi. Substituting these results particle caused by the fluctuating force during an 

in Eq. we get 



{mm} 



47r2 



dvif{vi) I d9ip{9i)sm(f> 
u 



4.2 The diffusion coefFicient 

If Av — Jjj* T{t')dt' denotes the increment of velocity of 
the test particle caused by the fluctuating force during £ 
interval of time i, we define the diffusion coefficient by 

((Z\w)2\ 



X sm 



■ sin(cjt) -t- (t)Cos{ijjt) 



m 

(33) This can be rewritten 



2t 



In the T — > limit, the spatial distribution of the particles 
^ can be approximated by 



D{t) 



1 

2i 



dt' / dt"{T{t')T{t")). 



(41) 



(42) 



p{9) = Mlo 



1/2 



Since the correlation function depends only on the time 
(34) interval \t" — we also have 



Using furthermore sin(a;) ~ a; in Eq. (|33|) for small 0, v and 
performing the Gaussian integrations, we finally obtain 



D{t) 



dt' / dt"{T{())T{t" ~t')). 



, , , , Affc2 v9 . , , MP 
{F{0)Fit)} ^ -^—sm{u;t) + 



47r2 



47r2 



1 



Setting T = t" — t' , we get 



cos(ajt). 
(35) 



1 /■* m' 



dT{T{0)T{T)), 



(43) 



(44) 
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or, equivalently, 



D{t) = - f dT f dt'{T{{))T{T)). (45) 



Finally, we obtain 



t 



D{t)^- / dT{t~T){mT{T)). 



(46) 



If the correlation function decreases sufficiently rapidly 
with time, taking t +oo, we obtain the Kubo formula 



D= {T{Q)T{T))dT. 
Jo 



(47) 



However, since in the present situation the temporal corre- 
lation function has an oscillatory behaviour, this formula 
is not applicable. According to Eq. (|^ . the diffusion co- 
efficient can be written D — Di — D2 where 



and 



i?i = / (^(0).F(r))dr, 

D.2 = \ j\mnr))rdT. 



Using Eq. we find that 
fc2 M 



D2 = 



P M 



47r2 



sin(a;i), 



47r2 /3w3 
Thus, we obtain 



smiwi) -I cosiut) 

^ ' LUt ^ ' Lot 



fc2 Ml- cos{ujt) 



(48) 



(49) 



(50) 



(51) 



(52) 



The diffusion coefficient is periodic and goes to zero at 
each time tn = {2tt /uo)n with n = 1,2, ... (see Fig.^. For 
t ^ 0, it behaves like 



D{t) 



fc2 M 



(53) 



On the other hand, D{t) for f ^ +00. 

4.3 The spatial correlation function 

Let us finally provide the exact expression of the spatial 
correlation function {T{9)J-{0')) in the case where the cor- 
relations between particles are neglected (as before). The 
effect of correlations is considered in |7l9j for the homo- 
geneous phase. The case of the inhomogeneous phase will 
be considered elsewhere. 




Fig. 1. The function D{t)/Do where Do ^ ''^ ^ 



Repeating the same steps as in Sec. 14.11 the spatial 
correlations of the fluctuating force can be written 



{mno')) = {F[e)F{e')} - -(Fie)){F{e')), (54) 



where 



{F{e)F{e')} 



47r2 



(55) 



and 



(Fie)} =B sine ^'-^hm^^o. 



27: loipB) 



(56) 



Using Eq. Q, we have 

X / sin(6l-6'i)sin(6l'-6li)e~''^™'^^M6'i. (57) 

Expanding the trigonometric functions and using the iden- 
tities 



/ sin6li cos6'ie-^^™"®id6li = 0, 

^0 



(58) 



/ sinH,e-P'''°'<'^d9i = ^h{pB), (59) 
Jo pB 

the integrals in Eq. (|57|l can be easily performed. Then, 
using Eqs. H56() and H54|l . we finally obtain 



47r2 



xIo{x) 



cos{e ~ e') 



1 - 



h{xf 2h{x) 
h{xY xh{x) 



sm a sm ( 



(60) 
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where we have set x = (3B. We note that, due to the 
inhomogeneity of the system, the correlation function of 
the fluctuating force is not a function of \9 — 0'\ alone. 

Let us consider particular cases. If we take 6' = 6, we 
see that {T{9Y) depends on 9 through a term sin^ 9. If we 
take 9 — Q and 9' — 6, we obtain 



PM Ii{x) 
Att^ xIq{x) 



COS( 



(61) 



For T >Tc (homogeneous phase) , we have a; = and Eq. 
reduces to 



(F(O)^^(^)) = 



"8^ 



■ COS( 



(62) 



We recover the result of I7'9' when the correlations be- 
tween particles are neglected. For T < Tc, using Eq. 
we find that 



kT 

(•^(0)^(0))-^ cos < 



(63) 



The dependence of this correlation function with the tem- 
perature is plotted in Fig. [21 Coming back to the function 
(|6U|) . and considering the homogeneous phase x = 0, we 
get 



I 

47r2 



9'). 



On the other hand, for T —> corresponding to x 
we find that 

kT 

{T{9)T{9')) = — cos(^ + 9') + 0{T^). 



kT 
2^ 



In particular. 

If we assume that 9 



{T') = f+0{T% 

2,71 

which coincides with Eq. (|4U|I when t — Q. 



cos(26i) + 0(t2). 
T as in Sec. EH then 



(64) 
-l-oo, 

(65) 
(66) 
(67) 



5 The friction force 

We shall now compute the frictional force 

(Ffr(i)) = P j dr J d9idviT{\ Q,t)T{Q -^l.t-r) 

xvi{t-T)p{9i)f{v,), (68) 

experienced by a test particle with prescribed trajectory 
given by Eq. (|16() . The expression (|68|l can be directly ob- 
tained from a linear response theory as done by Kandrup 
for the gravitational interaction. The friction arises 
as the response of the field particles to the perturbation 




Fig. 2. The normalized variance of the fluctuating force 
C(0) = if^iJ^iOf) at 61 = as a function of the tempera- 
ture when the correlations between particles are neglected. 



caused by the test particle as in a polarization process. 
Note that when the test particle is described by a distri- 
bution function P{9,v,t) [instead of having a prescribed 
trajectory] as in Sec^J the combination of terms H68() also 
enters in the friction term of equation ^ although the ex- 
pression is more complicated as it involves the history of 
the distribution function of the test particle. 

Since the fluctuating force averages to zero, Eq. (|68|l 
can also be written as 

(Ffr(<)) J dr J d9idviF{l ^ 0, t)T{0 ^l,t-T) 

xvi{t-T)p{9,)f{vi). (69) 

Taking the origin of times at t = 0, we shall first compute 
the quantity 



I{t) = J Fil ^ 0,0)F{1 ^ 0,t)vi{t) 

xp{9i)fivi)d9idvi, (70) 

which is related to the function !f'(i, r) defined in Eq. pifl . 
Using the fact that 



vi{t) = vi coa{Lut) — 9iuj ain{ujt), 
we can write I = Ii + I2 where 



(71) 



hit) = cos{ujt) I F(l -> 0, 0)F{1 -> 0, t)vi 

xpi9i)f{vi)d9idvu (72) 



I^it) = -Lusm{Lut) I F{1 -> 0,0)i^(l 0,t)9i 

xpiei)f{vi)d9idvi. (73) 
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The first integral can be rewritten 



h{t) ^cos(wi) J dvif{vi)vi J d9ip{9i)sm(j) 

— sin{ujt) + (j)cos{Lijt) . (74) 

CO J 

Using the same approximations as in Sec. 14.11 we obtain 
fc2 Me 



47r2 2w/3 



sin(2a;t). 



(75) 



Our asymptotic expansion is here valid to order T^/"^. The 
second integral can be rewritten 

hit) ^ ~^;^ujsva{ujt) j dvif{vi) j d9ip{9i)9i sincj) 



X sm 



— sin(a;t) + cos(a;t) , (76) 



leading to 

h{t) 



W- {cut) + "TT -5— sm(2wi) . (77) 
47r2 /3cj2 47r2 puj 



Summing these results, we get 
I{t) 



Since sin^(a;t) = [1 — cos(2a;t)]/2, the function I{t) is pe- 
riodic with pulsation 2lo, i.e. twice the pulsation of the 
orbiting particles. However, this function also contains a 
constant component which does not average to zero over 
a period. We have instead 



1 



m) = y / mdt 



87r2/3u;2 ' 



(79) 



where T here denotes the period. This implies that the 
effects of the friction accumulate with time. 

To obtain the complete friction force, we also have to 
calculate the quantity 



J(0- / Fil^O,0){Fil^O,t))v,{t) 
xpi0i)f(vi)deidvi. 



(80) 



Using 



(F(0^1,<)) 



k [vi 



2tt 



■ sm{Lut) + 61 cos{Lut) , (81) 



and performing the same approximations as before, we 
find that 

J{t) — f (j) (— sin ujt + 61 cos ujt^ 
47r^ J \ Lu ) 

X [v \ cos bjt — B\ijjs\mjjt) p{B\)f {y\)dQ\dv\. (82) 



20 



10 



-10 




10 

(Ot 



15 



20 



Fig. 3. The temporal behaviour of the components of the fric- 
tion force proportional to v and Q. 



By using the parity of the velocity and angular distribu- 
tions, we get 



J(<) = gsin(2c.t) I (J"^?-) 



p(0i)/(«i)d0id«i. 



(83) 



Finally, performing the Gaussian integrations, we find that 
J = 0. Therefore, the frictional force at time t is simply 
given by 



(Ffr(t)) = -/3 / /(-r)dr. 
Jo 



(84) 



Focusing on the component proportional to the velocity 
of the test particle (which contains the diverging contri- 
bution for t +00), we get 

(^ft-W) - -l^TJF I sin^(.^r)dr^~mv- (85) 

We note that, in the present situation, the friction coef- 
ficient is not given by an Einstein relation. Keeping only 
the diverging contribution coming from the non- vanishing 
averaged value of I{t) given by Eq. f7??| . we find that 



so that 



^, , k^M 1 kt 



{Fir{t)) = -^V{t). 



(86) 



(87) 



Therefore, the friction coefficient increases linearly with 
time. If we consider only that part (linear in t) of the 
friction coefficient, we find that it is related to the diffusion 
coefficient by 



2(1 — cos wt) ' 
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which resembles the Einstein relation. For i — s- 0, we find 
that ^{t) = D{t)(5. 

If we now consider the full expression of the friction 
force, we get 



sin2(jjt , - / „ X 
v\t + -{cos2ujt ~ 1) 

ZiU j 1 



(89) 



The temporal behaviours of the components of the friction 
force proportional to v and 9 are represented in Fig.|31 For 
t +00, we recover Eq. (|87|l and for t ^ 0, we have 



PM (v 



47r2 V3 2 



(90) 



Finally, with the quantities defined in Sees. ^ and |31 
the kinetic equation |(^ can be written 



dP , ,dP 
V— + (F)- 



dv 



dP 

'dt ' " 09 ' ' dv 
drl {T{t)T{t -t))—+ I3I{^t) 



dv 

■X.P{e{t - T),v{t - T),t - t). 



temporal correlations of Q{t). Furthermore, it is not clear 
whether Eq. Q conserves energy and whether it converges 
towards some form of equilibrium (Maxwellian or other) 
for large times. If we take for granted that the system 
must converge towards statistical equilibrium for t — > +oo, 
this implies that the kinetic theory may not be complete; 
one may have to relax certain simplifying approximations 
and consider next order terms in the expansion in 1/7V 
of the correlation functions. Alternatively, if we rely on 
the kinetic theory to give the justification (or not) of the 
statistical equilibrium state, one may conclude that the 
Boltzmannian distribution may not be reached (at least 
in a strict sense) for inhomogeneous systems. There might 
be dynamical anomalies preventing the system from reach- 
ing equilibrium. These are open questions left for future 
works. Note, however, that the preceding discussion does 
not favour other forms of entropy. The Boltzmann entropy 
remains the most fundamental even if the iJ-theorem may 
not be rigorously valid. 



6 Kinetic equation with time dependent 
friction 



(91) 6.1 General solution 



Using Eq. H40I) and keeping only the constant term in Eq, 
(CHI), we get 



dT 



dP 

2 cos(cjt) — 
av 



dP .^,dP 
^5^ + <^>a^ = 



Pvit)}P{eit 



The study of these kinetic equations is beyond the scope 
of the present paper. Let us simply state that the kinetic 
theory of inhomogeneous systems with long-range interac- 
tions is far from being completely understood. In partic- 
ular, it is not clear whether Eq. H92|l . and more generally 
Eq. {Tj), satisfies an i7-theorem. Indeed, starting from the 
general kinetic equation Q, introducing the Boltzmann 
entropy S = — J / In fdOdv and using standard methods, 
we can put the rate of entropy production in the form 



The kinetic equation H92|l is complicated due to non Marko- 
vian effects and spatial inhomogeneity encapsulated in the 
2 advection term. In addition, one striking novel aspect of 

k M 1 d our study is the time dependence of the diffusion coefR- 

Stt^ /Jw^ dv cient and friction force. Therefore, in a first attempt to 
investigate the effects produced by such terms, we shall 
r),w(t — r),i — r). consider a simpler kinetic equation where we neglect non- 
Markovian effects and advective terms altogether but keep 
(92)the time dependence of the diffusion and friction coeffi- 
cients. In that case, Eq. 192|l reduces to 



dP d f dP 



(95) 



2 



dOdvdOidvi — — 
fji 



dTQ{t)G{t,t^T)Q{t 



(93) 



where 



where D{t) and ^{t) are given by Eqs. and More 
generally, we shall consider equations of the form (|95() for 
arbitrary functions of time D{t) and S,(t). An interesting 
aspect of the problem is that such equations can be solved 
analytically. We note that when D and ^ are constant, we 
get the familiar Kramers equation. Its stationary solution 
is the Maxwellian P^ = Ae~^^ 1"^ provided that D and ^ 
are related to each other by the Einstein relation ^ = D(i. 
We shall now consider the case where D and ^ depend on 
time. By redefining time such that dt' = D{t)dt, we are 
led to consider, without loss of generality, equations of the 
form 



Q{t) 



T{l^Q,t)—+T{Q^l,t) — 

OV OVi 

xf{e,v,t)f{0,,v,,t) 



(94) 



dP d f dP , , 



(96) 



In the course of the calculations, we have inverted the 
dummy variables 9,v and 9i,vi and taken the half sum 
of the resulting expressions (see, e.g., ^|). We see that 
the iJ-theorem S* > is not granted. This depends on the 



where D and 7 are constant parameters. 

Taking the Fourier transform of Eq. (|96|l with the con- 
ventions 



P{v) = / PiOe'^-^di, F(0 = / P{v)e 



dv 
2^' 
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(97) g{v) = G{v/2x{t)) 



and using the relation 

P{v)v = j vP{()e'^'"d^ 

4/p(0|(e^^"K = 4/fe^^"^e, (98) 



we get 



(99) 



We introduce the change of variables 

f{y,t) = P{Hi{t)y,t), e = i/i(% (100) 
and choose the function Hi{t) such that 



Substituting Eq. ((TUUI) in Eq. (jHS, we find that f{y,t) 
satisfies 



Let H{t) be the primitive of h{t) such that 



H{t) = / h{T)dT. 



(102) 



(103) 



Then, we choose Hi, solution of Eq. pUl|) . such that 



Hi{t) 



(104) 



By convention, 7J(0) = and i?i(0) = 1. Equation ((TU^ 
can be integrated leading to 

/(2/,^) = /(2/,0)e-^^W^^ (105) 
where we have defined 

H2{t)^D [ HiirfdT. (106) 
Jo 

Returning to original variables, we obtain 



Hi{t) 



(107) 



where we have defined 



Hi{tf JO 



D e''<Wr)-mt)]dr. (108) 



Defining 



(109) 



m = 2x{t)G{2xm 



where 



G{z) = e- 



G(0 



20F 



we can rewrite Eq. H1Q7|I in the form 



(110) 



(111) 



(112) 



Taking the inverse Fourier transform, we can express the 
solution of Eq. H96|l as a convolution 

P(^;,^) = ^A^ j q{v - v')g{v')^, (113) 



(101) or, equivalently 



Piv,t) = ^ e--'Po[Hi{t){v-2x{t)x)]dx. 



(114) 

By direct substitution, we can check that Eq. (|114|l is in- 
deed solution of Eq. 1)96(1 . For 7 = 6 (pure diffusion) we 
find that 



P{v,t) 



er"" Po{v-2x^/Dt)dx. (115) 



For I? = 0, Eq. (|114|) reduces to 

Piv,t) = Hi{t)Po[Hi{t)v] 



(116) 



This can also be obtained by noticing that for D ~ 0, Eq. 
(|96|l is an equation of continuity. The equation of charac- 
teristic is dv/dt = —^h{t)v which can be integrated into 
v{t) = voe-^"^*\ Writing P(v,t)dv = PoMdvo, we fi- 
nally get Eq. l(TTC|l . 

If ^'o(''^) = vi^~^o) is a step function with rj{v—Vf)) — 1 
for V < vo and rj{v ~ vq) — for v > va, we find that 



P{v,t) = Hi{t)$ 



V -vo 



/Him 



2x(i) y 



where 



1 



+ CXD 



e ^ dy. 



(117) 



(118) 



Alternatively, if Poiv) — 8(v — vq), we obtain 

1 (^-^,o/gi(t))^ 

■ ——7% 



P{v,t) = 



(119) 



P.H. Chavanis: Strongly inhomogeneous HMF model 



11 



6.2 The case h{t) = 1 

In the ordinary situation where the friction coefficient is 
constant {h{t) = 1), we get 



Hit)=t, Hi{t) = e 



7* 



(120) 



Equation pi7(l then takes the form 



(121) 



(122) 



which behaves like 



Piv,t) ~ e'^*'l> 



V — VqC 



-ft 



(123) 



for t +00. On the other hand, Eq. (|119|l takes the form 



1 



■ T(^-^oe-V')" 

e 2I,(l-e-^7t) ^ (124) 



27r7:>(l - e-27t) 
which tends to the Maxwelhan for t ^ +oo. 

6.3 The case h{t) = t 

In the case where the friction coefficient increases hnearly 
with time {h{t) = t), we get 



(125) 



H2it)=D e'^^'dT, xHt)=De-'''l e^'^'dT. 
Jo Jo 

(126) 

The last function can be written x^(i) = 'y=4'iy/lt) where 



^(t) = e-' / da:, 
Jq 



(127) 



is the Dawson integral. It behaves like 4>{t) ^ l/(2t) for 
t +00 and like (j){t) ~ t for t 0. Equation pi7|) then 
takes the form 



P(.,t) = e^*V.^|^^^|Ll2l,, 



(128) 



which behaves like 



/ 7t(« - «oe-^*'/2) 



/2D 



(129) 




Fig. 4. The function P{v,t) with 7 = D = 1 and vg — i for 
t = 0.01, 0.1, 1, 1.5 starting from a step function ri{v — vq). 



for t —* +00. The time evolution of the profile P{v,t) is 
represented in Fig. 01 On the other hand, Eq. (|119|l takes 
the form 



P{v,t) = 



1 



1/2 



For t — > +00, it behaves like 
/ + \ 1/2 



(130) 



P(«,t) 



V27ri:> 



7* / -Tt^/2\2 



, (131) 



and eventually tends to 5{v). The system collapses to a 
Dirac in velocity space due to the divergence of the friction 
coefficient for t +00. Alternatively, for i ^ 0, the initial 
distribution 5{v — vo) deforms itself according to 



P{v,t) 



1 („-„oe-"'*V2)2 

=e -i-D* 



(132) 



The time evolution of the profile P(w, i) is represented in 
Fig. El 



6.4 Another example 

For completeness, we also provide the analytical solution 
of the equation 



dP d fdP 



(133) 



where D and 7 are constant parameters. Writing this 
equation in the form 



dP , , ,dP d'^P 



(134) 
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Fig. 5. The function P{v,t) with ^ — D = 1 and wo = 4 for 
t = 0.05, 0.1, 0.2, 0.5, 1, 2, 3, 5, 10 starting from 5{v - vo). 



we see that the second term is siniilar to an advection 
(in velocity space) by an effective velocity field V{v, t) — 
—^h{t). Let Vf{t) be a solution of dvf/dt = —jh{t). We 
take Vf{t) = —jH{t) where H{t) denotes a primitive of 
h(t) with H(0) — 0. Then, we define z — v — Vf{t) and 
P{v, t) = (j){z, t). Substituting these relations in Eq. 1134(1 . 
we find that (t>{z,t) satisfies the diffusion equation 



dt 



D 



92 

dz 



(135) 



Using for example the results of Sec. 16.11 the general so- 
lution of this equation is 



iz,t) 



1 



e'"" Mz~'2vDtx)dx. (136) 



Returning to original variables, we get 
2 



Poiv + jH{t) - 2VDtx)dx, 

(137) 



where Po(''^) = P{v,0) is the initial value of the probability 
distribution. If Po(v) = 7]{v — vq) is a step function, we 
find that 



P{v,t) = ^ 



lH{t) \ 



2VDt J 
Alternatively, if Po{v) = 5{v — wq), we obtain 



P{v,t) 



7 Conclusion 



1 



v'47rDi 



e ■lo' 



(138) 



(139) 



We have presented first elements of kinetic theory appro- 
priate to inhomogeneous systems with long-range inter- 
actions. Explicit results have been obtained for the HMF 



model in the low temperature limit T ^ where the par- 
ticles perform harmonic motion around = with the 
same frequency. These results show that the description 
of strongly inhomogeneous systems is very different from 
what we are used to in ordinary kinetic theory. In particu- 
lar, the Kubo formula is not valid in its usual form because 
the auto-correlation function of the force does not decrease 
sufficiently rapidly for large times. On the contrary, it has 
a striking oscillatory behaviour that averages to zero over 
a period. On the other hand, the Einstein relation is bro- 
ken. Indeed, the friction coefficient depends on time and 
diverges linearly as t —t +00 due to the cumulative effects 
of friction as the particle undergoes several oscillations. 

These results are strikingly different from those ob- 
tained in the homogeneous phase T > Tc oi the HMF 
model where the particles follow linear trajectories with 
constant velocity in a first approximation. In that case, the 
Kubo formula and the Einstein relation are recovered 
It would be interesting to consider the case of the weakly 
inhomogeneous HMF model. This could be studied pertur- 
batively by expanding the relations H13(l - p4(l close to the 
critical point T — > where B ^ Q. Similar expansions 
have been previously considered in '7' in other circum- 
stances. It will be necessary to distinguish between open 
trajectories that are similar to the situation for T > Tc 
from closed trajectories that are similar to the situation 
for T ^ 0. It would also be interesting to investigate how 
collective effects can modify these results (as in the case 
T > Tc). In particular, close to the critical point the corre- 
lation time diverges which complicates the kinetic theory. 

We stress that our approach does not provide a full 
kinetic theory of the inhomogeneous HMF model. Our 
main motivation was to mention the difficulties with such 
a kinetic approach due to Markovian effects and spatial 
inhomogeneities and to show the differences with more 
familiar kinetic theories. Even if we do not obtain many 
explicit predictions, an interesting aspect of our study is 
to show that the usual kinetic relations (Kubo formula, 
Einstein relation,...) can break down for strongly inho- 
mogenous systems with long-range interactions such as 
the HMF model. Finally, we would like to mention that it 
is important to develop a kinetic theory for the inhomoge- 
neous HMF model because there are situations of physical 
interest where the system is spatially inhomogeneous. On 
the other hand, the bath can have a complicated phase- 
space structure with a hierarchical cluster size distribution 
as exemplified in the numerical simulations of Rapisarda 
& Pluchino 121 ■ The kinetic theory should take into ac- 
count these effects. 
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